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Bivariate extreme-value distributions have been used in modeling extremes in environmental 
sciences and risk management. An important issue is estimating the dependence function, such 
as the Pickands dependence function. Some estimators for the Pickands dependence function 
have been studied by assuming that the marginals are known. Recently, Genest and Segers [Ann. 
Statist. 37 (2009) 2990-3022] derived the asymptotic distributions of those proposed estimators 
with marginal distributions replaced by the empirical distributions. In this article, we propose 
a class of weighted estimators including those of Genest and Segers (2009) as special cases. We 
propose a jackknife empirical likelihood method for constructing confidence intervals for the 
Pickands dependence function, which avoids estimating the complicated asymptotic variance. 
A simulation study demonstrates the effectiveness of our proposed jackknife empirical likelihood 
method. 
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1. Introduction 

Let (Xn,Xi2), ■ ■ ■ , (X n i,X n 2) be independent random pairs with common distribution 
function F and continuous marginal distributions Fi(x) = F(x, 00) and F 2 {y) = F(oo, y). 
Then the copula of F is defined as 

C(x,y)=P(F 1 (X n ) < x,F 2 (X 12 ) < y). 

When C^u 1 /*, v 1 ^) = C{u,v) holds for all u,v € [0, 1] and t > 0, C is called an extreme 
value copula and is determined by the Pickands dependence function, A, through the 
equation 



C(.,,)=exp{log M A(^)} 



(1.1) 
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for all (u, v) <E (0, l] 2 \ {(1, 1)}, where A is a convex function and satisfies max(t, l — t)< 
A(t) < 1 for all < t < 1 (see Pickands [16] and Falk and Reiss [6]). 
Write Y i:j = - U>»{ F,\ X ,/>} for i = 1, . . . , n, j = 1,2 and 



H n { z ) = -Yl( Ya <z). 



We denote u A v = imii(u,v) and u V v = max(u,u) throughout. Estimators for the 
Pickands dependence function A(t) when the marginal distributions Fj,j = 1,2 arc known 
have been proposed by Pickands [16], Deheuvels [5], Hall and Tajvidi [10], and Caperaa, 
Fougeres and Genest [3], defined as 

A P (t) 
A D {t) 
A HT (t) 



£1LiOV*}aOV(i-*)} : 



EIUOV*} a {Y a /(i t)} tj2ti Yn - (i - 1) ELi *fe 

EUP^iWES Y n )} A {r,y i2 /((l - t) J2U Y 2)} ' 



A^(t) = e«p{A(t) /*M!^£d*-(l-A(t)) f^T^dz 

respectively, where A(t) € [0,1] is a weight function and ^4 P (i) and .^(i) are cor- 
responding limits when t = or 1. When the marginal distributions are unknown, 
similar nonparamctric estimators can be obtained by replacing the marginal distri- 
bution Fj by the corresponding empirical distribution F n j(x) = ^ E"=i ^(^ij — x ) or 
F nj (x) = ^j-ELi ^ x )- Wc denote these estimators as A p {t), A D (t), A HT (t) and 

A CFG (t). Recently, Genest and Sogers [8] showed that A p (t), A D {t) and A HT (t) have 
the same asymptotic distribution as 

A P (t) 



and that A CFG (t) with A(i) = t has the same asymptotic distribution as 
A CFG (t) = exp j-7 - i - *)) A | , 

where 7 = — J °° log(x)e _2: da; is the Eulcr constant and 

Zij = -log{F nj (Xij)} for i = l,...,n,j = l,2. 

Moreover, Genest and Segers [8] derived the asymptotic distributions of A p (t) and 
A (t) by noting the following important relationship: 



1 

-1/ 



A F (t) = { / M - 1 C„(u 1 - t ,u t )dw 
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and 

A CFG (t) = + ^{Cniu 1 -*, «*) - 7(u > e- 1 )}^^)}- 1 du) , 

where 

1 " 

<?„(«,«) =Tl(F„i(Iii) < u,F n2 {X l2 ) < v). 

i=l 

In this article, we propose a class of weighted estimators including A p (t) and A CFG (t) 
as special cases. We provide details in Section 2. In Section 3 we propose a jackknife 
empirical likelihood method to construct confidence intervals for the Pickands dependence 
function. Unlike the normal approximation method, this new method does not need to 
estimate any additional quantities, such as asymptotic variance. In Section 4 we report 
a simulation study conducted to examine the finite sample behavior of the proposed 
jackknife empirical likelihood method. We provide proofs in Section 5. 



2. Weighted estimation 

It follows from (1.1) that 

C{u 1 - t ,u t )=u A{ - t) for all u€ [0,1] and alliG [0,1], (2.1) 

which motivates the estimation of A(t) by minimizing the following weighted distance 
with respect to a > 0: 

f {C' n (u 1 " t ,it t )-M Q } 2 A(it,f)du, 
Jo 

where A(tt, t) > is a weight function. Under some regularity conditions, the foregoing 
estimator is the solution of a to the equation 

/ {C*„(w 1- *, u*) - u a }u a {- log(u)}A(u, t) du = 
Jo 

for a > 0. This is a special case of the proposed M-estimators and Z-estimators of Bucher, 
Dette and Volgushev [2]. Noting that u a (—\ogu)A(u,t) = C(u 1 ~ t ,u t )(—\ogu)A(u,t) and 
X(u,t) is any weight function, we propose treating (7(u 1 ~ t ,u*)(— logu)A(u,t) as a new 
weight function. This leads us to estimate A(t) by solving the following equation with 
respect to a > 0: 

/ {C , Il (u 1 -*,w t )-M Q }A(u,i)dM = 0, (2.2) 
Jo 

where X(u,t) > is a new weight function. We denote this new estimator by A™(t;X). 
When X(u,t) is taken as u" 1 or {-iilog(u)}" 1 , A%(t;\) becomes A p (t) or A CFG (t). 
Thus, the foregoing class of estimators includes the known estimators in the literature 
as special cases. 
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Write g(a) = f Q {C n {u 1 ~ t , u t ) — u a }X(u,t) du. Because u a is a decreasing function of 
a for each fixed u £ [0, 1], g(a) is an increasing function of a for each fixed t. Moreover, 
tjf(O) < and g(oo) > when n is sufficiently large. Thus, (2.2) has a unique solution 
A™(t; A) for each large n and t £ [0, 1]. Note that this unique solution might not satisfy 
that max(t,l-t)<A%(t;X)<l and i™(0; A) =i™(l; A) = 1. 

Let W(u, v) denote a tight Gaussian process with mean 0, covariancc 

E{W(u 1 ,v 1 )W(u 2 ,v 2 )} = C(ui Au 2 ,V! Av 2 ) - C(ui,vi)C(u 2 ,v 2 ), 

and W(u, 0) = W(0, v) = W(l, 1) = for all u,v£ [0, 1]. The asymptotic distribution for 
the proposed estimator A™ (i; A) is given in the following theorem. 

Theorem 2.1. Suppose that ^;C(u,!)), -^C{u,v) and Q^ dv C(u,v) are defined and 
continuous on the sets J-\ = {(u,v): < u < 1 and < v < 1} , T-z = {(u,v): < u < 
1,0 < v < 1} and J-3 = {(u,v): < u < 1,0 < v < 1}, respectively. Also assume that for 
each fixed t G [0, 1] the function X(u,t) > is continuous and not equal to as a function 
of u G (0,1). Furthermore, assume that 



d 2 



dv 2 

d 2 



< 



M 



du dv 



C(u,v) 
C(u,v) 



< 



u(l — u) 

M 
v(l — v) 
M 

< 



M 



u(l — u) v(l — v) 



for (u,v) € T\, 
for (u,v) e T 2 , 
for (u, v) e 73 



for some constant M > 0, A'(t) is continuous on [0, 1], and there exist 6± > and 6 2 € 
[0,1/2) such that 



sup y/ri I 

0<i<l Jo 



-l/((l-t)Vt) 



u 1/2 X{u,t)du^0, 



sup \fn \ (1 — u)X(u,t) du — > 0, 

0<t<l J(n/(n+l)) 1 /« 1 - t > vt > 

(n/(n+l)) 1/((1 ^ t)Vt) 

X(u,t)du^0, 

( rl+ l)-i/((i-t)vt) 



sup n 

0<t<l 



-l/4+5i 



sup 

0<i<l Jo 



sup 

0<t<l Jo 



sup 

0<t<l Jo 



sup 



{«(i-*)v*( 1 _ u (i-*) vt )}^A( u ,t)d M < 00, 
u (i- t) vt-(i- t ) u (i^ 2(1 _ u i-t)& A(u>t ) du < ^ 
u (i-*)v*-t u t& _ u * j« a A ( Uj i) du < 00; 



(2.3) 



L 0<t<l Jl/2 



(— log«)A(u, t)dt < 00. 
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Then as n — > oo, swp 0<t<1 \A™(t] \) ~ A(t)\ =o p (l). Moreover, suppose that X(u,t) is 
continuous in (0, 1) x [0, 1] and 

[X(u,h) - X(u,t 2 )[ < \h-t 2 \ So X Q (u), ti,i 2 6[0,l],wG(0,l) 

for some constant 5q > and function Xq(u),u € (0, 1), where Xq(u) satisfies that 

•1/2 



u a Xo(u) du < oo, 



(1 — u a )Xo(u) du < oo 



1/2 



for all a > 0. Then, asn->oo, y / n{A^'(t; A) — converges to B(t) in C ([0,1]), where 



Bit) 



C(M 1 -',u')A(u,t)log(u)dM 



x / {^(u 1 -*,^) - C^m 1 -*,^)^^ 1 "*, 1) - C 2 (u 1 - t ,u t )W(l,u t )}X(u,t)du, 
Jo 

d(u,v) = J^C(u,v) and C 2 (u,v) = J^C(u,v). 

Remark 2.1. Theorem 2.1 still holds when condition (2.3) is replaced by 



sup \/n u x / 2 X(u,t) du — >• 0, 

o<t<i Jo 

(1 - u)X(u,t) du -> 0, 



sup 


sjn 


0<t<l 




sup 


\fn 


0<t<l 




sup 


n- 1 


0<i<l 




sup 


f 


k 0<t<l 


Jo 



'(n/(n+l)) 2 

,(n/(«+l)) 2 

I ' i+Sl / A(u,f)du->-0, 

J(n+1)- 2 

u S2/2 (l-u) 52 X(u,t)du<oo 



for some 5i > and 5? E [0, 1/2). This follows from the proof of Theorem 2.1 by replac- 
m s Jo ' J(n+i)- 1 /C( 1 -')vt) ana J( n /(„ +1 ))i/((i-t)vt) m 1,0. oj Dy j Q , 

j{n+i n )-* )] and /(l/Cn+i))^ respectively. 

Remark 2.2. A common approach to choosing X(u,t) is to minimize the asymptotic 
variance of A™(t;X). This is difficult to do analytically. Linear combinations of some 
known estimators can be considered instead. For example, suppose that the weight func- 
tions Xi(u), . . . , X q (u) give the corresponding estimators A™ 1 (t), . . . , A™ q (t). Define the 
class of new weight functions as 



Jb = I X(u, t): X(u,t) = J2 ai (t)Xi (u) , Oi (t) > 0, . . . , a q (t) > 0, ^ a t (t) = 1 
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Then one can choose a[s to minimize the asymptotic variance of A™(t;\) in this class 
.7-0, which results in explicit formulas for a^s. 

An example. Assume that X(u,t) = \ogu)~ q ^ for some q(t) £ [0, 1]. Then A p (t) 

and A CFG {t) correspond to q(t) = and q(t) = 1, respectively. When q(t) < 1, we can 
write 

/•i 

{C n (u 1 - t ,u t ) -u e }X(u,t)du 

= ~Z —TTV / {Cntf-W-u'M-togu) 1 - 9 ® 

1 ~ K 1 ) JO 



1 



i - q(t) Jo 



1 r 1 fjiW- 1 r°° 

- {-loguf-^AC^-'M)-- u^e-Mi 

1 — QK 1 ) Jo 1 iv-) Jo 



_1 I lyj^L ^2 1 1 _ r (t)-i r(2 _ 



where the ZLs are as defined in Section 1. Thus, 



i-(t; A) = exp<j - ( log( i £ A ^] ' j - lo g r(2 - g(i)) ) /(I - 9 (t)) 



for < q(t) < 1. Note that when q(t) = 1, the foregoing expression is defined as the 
limit, which becomes the same as A CFG (t). In particular, we propose to choose q(t) = 
mm{A CFG (t), 1} and denote the resulting estimator by A™(t). To compare this new 
estimator with A c (t), we draw 1000 random samples with size n= 100,1000,5000 
from a Gumbel copula with A(t) = {t 9 + (1 — t) 6 } 1 / 6 , a Hiisler-Reiss copula with A(t) = 
(1 - t)$(6 + ± log ±=£) + 1$> {9 + ± log j^), and a Tawn copula with A(t) = 1 - 6t + 9t 2 , 
where &(x) denotes the distribution function of N(0, 1). Figure 1 plots the ratios of the 
mean squared error of A™(t) to the mean squared error of A CFG {t) for t = 0.1, 0.2, . . . , 0.9, 
and shows that the new estimator has a smaller mean squared error than A GFG (t) in all 
of the cases considered. 



3. Jackknife empirical likelihood method 

In this section, we consider interval estimation for the Pickands dependence function 
A(t), which plays an important role in risk management since one may be concerned 
with interval estimation for C(u, v) at some particular values of u and v. Note that an 
interval for A(t) can be easily transformed to an interval for a monotone function of A(t). 
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Husler-Reiss copula with theta=0.5 



Husler-Reiss copula with theta=2 



n=100 

n=1000 

n=5000 




0.4 0.6 
t 

Tawn copula with theta=0.25 



Tawn copula with theta=0.75 



n=100 

n=1000 

n=5000 















n=100 






n=1000 






n=5000 





Figure 1. Ratios of the mean squared error of the new estimator A™{i) to that of A CFG {i) for 
t = 0.1, 0.2,..., 0.9. 



Moreover, these two intervals have the same coverage probability, but different interval 
lengths. Because the upper tail dependence coefficient can be written as a monotone 
function of ^4(1/2), an interval can be constructed via an interval for A(l/2). 

An obvious approach to constructing an interval for A{t) is to use the normal approx- 
imation method based on any one of the estimators for A(t). Because the asymptotic 
distribution of any one of the estimators for A(t) depends on its derivative A'(t), the nor- 
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mal approximation method requires estimating A'(t) first. In an alternative approach to 
constructing confidence intervals, the empirical likelihood method has been extended and 
applied in various fields since Owen [13, 14] introduced it for construction of a confidence 
interval/region for a mean. (See Owen [15] for an overview.) An important feature of the 
empirical likelihood method is its property of sclf-studcntization, which avoids estimat- 
ing the asymptotic variance explicitly A general approach to formulating the empirical 
likelihood function is based on estimating equations, as in Qin and Lawless [17]. 

Because our proposed weighted estimator is defined as the solution to equation (2.2), 
the method of Qin and Lawless [17] may be applied directly by defining the empirical 
likelihood function as 

{n n 
JJ(npi): pi >0,...,p n > 0,y^ = 1, 
»=1 i=l 
n .\ 

Y^Pi / {I(F nl (X ll )<u 1 - t ,F n2 (X l2 )<u t )-u 9 }X(u,t)du = 
i=i J ° 

However, this method cannot catch the variation introduced by the marginal empiri- 
cal distributions. In other words, the limit is no longer a chi-squared distribution. In 
general, the nonlinear functional must be linearized by introducing some link variables 
before the profile empirical likelihood method is used. (See Chen, Peng and Zhao [4] for 
details on applying the profile empirical likelihood method to copulas.) Unfortunately, 
this linearization idea is not applicable to the estimation of A(t). Recently, Jing, Yuan 
and Zhou [11] proposed a so-called "jackknifc empirical likelihood" method to construct 
confidence intervals for U-statistics. More specifically, these authors proposed applying 
the empirical likelihood method to jackknife samples, which could result in a chi-squared 
limit. Motivated by Gong, Peng and Qi [9]'s study of the use of a smoothed jackknifc 
empirical likelihood method to construct a confidence interval for a receiver operating 
characteristic curve, one needs to work with a smoothed version of the left-hand side of 
(2.2). The reason for smoothing is to separate marginals from the copula estimator when 
expanding the jackknife empirical likelihood ratio. In this work, we used the smoothed 
empirical copula of Fcrmanian, Radulovic and Wegkamp [7] , defined as 

C-(u>-, „•) = I ±K (-^WO) K (-*ff<*°>) , 

where K(x) — J_ k(s) ds, k is a symmetric density function with support [—1,1], and 
h = h{n) > is a bandwidth. Based on this smoothed estimation, a jackknifc empirical 
likelihood function can be constructed as follows. Put F n j t -i(x) = ^ Y^i=\ i^i^i-^-lj — x ) 
for j = 1,2 and i = 1, . . . ,n, 
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for i= 1, . . . , n, and define the jackknife sample as 

Vi(u,t) =nC s n (u 1 - t ,u t ) - (n- lJC'^u'-V) 

for i = 1, ... ,n. 

We next apply the empirical likelihood method based on estimating equations of Qin 
and Lawless [17] to the foregoing jackknife sample. This gives the jackknife empirical 
likelihood function for 9 = A(t) as 



L{6) = sup< JJ(ripi): pi > 0, . . .,p n > 0,^p, = 1, 



.i=i 



i=l 



l-b„ 



{V t (u, t) - u e }X(u, t) du = }, 



where a n > and b„ > 0. Note that we use J 1 6 " instead of J Q in defining the foregoing 
jackknife empirical likelihood function, to control the bias term and allow the possibility 
of A(0, t) = oo and A(l, t) = oo. 

By the standard Lagrange multiplier technique, we obtain the log jackknife empirical 
likelihood ratio as 

n 

1(9) = -21ogL(0) = 2^1og{l + /3Q i (0)}, 



where 



and f3 = (3(9) satisfies 



Qi(0) 



1-6- 



- u e }A(u,*) dw 



1 ™ 

71 ^ 1 

»=1 



Qi(0) 



■PQi(0) 



(3-1) 



Theorem 3.1. Suppose that -J^2C(u,v) , J^j2C(u,v) and q^q v C(u, v) are defined and 
continuous on the set J-3 = {(u, v), < u < 1 and < v < 1} and 



y 2 



9it dv 



C(u,v) 
C(u,v) 



< 



M 



u(l — u) ' 



M 

< — A 



dv 2 
M 



C(u,v) 



< 



M 



v(l — v) 



u(l — u) v(l — v) 



for (u,v) € J-3 and some constant M > 0. Let t denote a fixed point in (0, 1). Assume that 
the function X(u,t) > is continuous and not identical to as a function of «£ (0, 1), 
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A'(s) is continuous at s = t, and 

h = h(n) — > 0, nh — >oo, 
a n 0, b n -> 0, /i/ a„ ->• 0, /i/6„ 

„-l/4+<5i 



0, 



l-6„ 



A(w, i) du — > /or some <5i > 0, 



w 52 (l-w' 52 )A('u,t)dw<oo for some 6 2 E [0,1/2), 



\fnh 



i-6 n 

2 / „.-3/2 



ti" ;i/2 A(u,t)dw^0, 



(3.2) 



i-&„ 



V^/i 2 / {logw} ' 

• 1-b 



1 



\fnh J a 



rl-b„ 
,-3/2 / „ -2 



~ 1 w~ 3/2 A(u,t)du^0, 
u~ 1 X(u,t) du^O, 
u~ 2 \(u,t) du ->• 

as n — > oo. TTien Z(Ao(t)) A x 2 (l) as oo, where Ao(t) denotes the true value of A(t). 

For any fixed f € (0, 1), based on the foregoing theorem, a jackknife empirical likelihood 
confidence interval for Ao(t) with level 70 can be constructed as 

/ To (t) = {0: l(0)< X * o }, 
where x 2 is the 70 quantile of X 2 (l)> as follows: 

Remark 3.1. (i) When X(u,t) = {— ulogu} -1 , we have sup 0<u<1 X(u, t) = 00. We can 
choose 



-1/3 



a n = d±n a , b„ = d 2 n , h = d^n 

for some di, 0^2,03 > 0, < a < 1/9, and < b < 1/6. 

(ii) When sup 0<u<1 \(u,t) < 00, we can choose 

a n =d\n~ a 1 b n =d2n~ h 1 h — d^n^ 1 ^ 3 

for some di, d 2 ,ds > 0,6 > 0, and < a < 1/3. Here we fix the rate for h because the 
optimal rate for the bandwidth in smoothing distribution estimation is n~ 1//3 . 

(iii) Theorem 3.1 still holds when a n — ^ a e (0, 1/2) and b n — > b € (0, 1/2] as n — >• 00. 



4. Simulation study 

In this section we examine the finite-sample behavior of the proposed jackknife empir- 
ical likelihood method based on X(u,t) = \ogu)~ min { A in terms of cov- 
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Table 1. Empirical coverage probabilities are reported for the proposed jackknife empirical 
likelihood confidence interval (JELCI) based on X(u,t) = logu)~ mm{A OFa (t),i} ^ anc j tne 

confidence interval based on the multiplier method for A CFG {t) (MCI) with nominal levels 0.9 
and 0.95 





Level 0.9 


Level 0.9 


Level 0.95 


Level 0.95 


(n, t, Copula, 9) 


JELCI 


MCI 


JELCI 


MCI 


(100,0.1,Gumbel,2) 


0.604 


0.276 


0.639 


0.366 


(100, 0.1, Hiisler-Reiss, 0.5) 


0.845 


0.566 


0.899 


0.655 


(100, 0.1, Tawn, 0.25) 


0.817 


0.571 


0.872 


0.670 


(100,0.5,Gumbel,2) 


0.871 


0.722 


0.941 


0.784 


(100, 0.5, Hiisler-Reiss, 0.5) 


0.888 


0.715 


0.941 


0.802 


(100, 0.5, Tawn, 0.25) 


0.886 


0.750 


0.941 


0.825 


(100,0.8,Gumbel,2) 


0.841 


0.531 


0.889 


0.599 


(100, 0.8, Hiisler-Reiss, 0.5) 


0.889 


0.646 


0.947 


0.758 


(100, 0.8, Tawn, 0.25) 


0.884 


0.677 


0.938 


0.758 


(1000,0.1,Gumbel,2) 


0.888 


0.655 


0.935 


0.740 


(1000, 0.1, Hiisler-Reiss, 0.5) 


0.892 


0.813 


0.942 


0.883 


(1000, 0.1, Tawn, 0.25) 


0.900 


0.820 


0.957 


0.891 



erage probability and compare it with the method based on the asymptotic distri- 
bution of A GFG (t). For computing the coverage probability of the proposed jackknife 
empirical likelihood method, we choose k(x) = j|(l — s 2 ) 2 /(|a;| < 1), h = 0.5n -1 / 3 , 

a n = b n = 0.1, X(u,t) = u _1 (-logw) _min ^ CFG(t) ' 1} and use the R package "emplik" 
(see Zhou [19]). To compute the confidence interval based on the asymptotic distri- 
bution of A CFG (t), we use the multiplier method proposed by Kojadinovic and Yan [12]. 
More specifically, we use cq. (7) of Kojadinovic and Yan [12], with TV = 500 and 
{Z^: i = 1, . . .,n, k = 1, . . ., N} as independent random variables from iV(0, 1), to cal- 
culate the critical points of the asymptotic distribution of ^/n{A CFG (t) — A(t)}. We 
do not use a larger N, because this multiplier method is computationally intensive. 
A comparison study on bootstrap approximations has been reported by Biicher and 
Dette [1]. 

We draw 1000 random samples with size n = 100, 1000 from the Gumbel cop- 
ula, the Hiisler-Reiss copula, and the Tawn copula with Pickands dependence func- 
tions specified at the end of Section 2. Table 1 reports the coverage probabili- 
ties at levels 0.9 and 0.95 for t = 0.1, 0.5, 0.8. These show that (i) the proposed 
jackknife empirical likelihood method gives much more accurate coverage probabil- 
ities than the multiplier method based on the asymptotic distribution of A G (t), 
and (ii) our proposed jackknife empirical likelihood method performs poorly for the 
boundary case t = 0.1 when n = 100, but its performance improves as n becomes 
large. 
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Proof of Theorem 2.1. Define 



a n (u, v) = v^j i £ I(F 1 (X u ) < u, F 2 (X i2 ) < v) - G{u, v) j . 

Then, from Proposition 4.2 of Segcrs [18] and Theorem G.l of Genest and Sogers [8], it 
follows that 

sup \y/n{C n (u,v) - C(u,v)} - a n (u,v) + Ci(u,v)a n (u, 1) + C 2 (u,v)a„(l,v)\ 

0<u,v<l 



: 0(n- 1/4 (logn) 1/2 (loglogn)) 



and 



a n (u,v) d W(u,v) 



(uAv) s (l -uAv) 5 (u A v) 5 (l - u A v) 6 

in the space /°°([0, l] 2 ) of bounded, real- valued functions on [0, l] 2 for any 5 £ [0, 1/2), 
where W is defined before Theorem 2.1. By the Skorohod construction, there exists a 
probability space carrying C*,a^,W* such that 



(C*,a* n ) = (C n ,a n ), W* = W, 
sup \y/n{C*(u,v) - C{u,v)} - a* n {u,v) + d(u,v)a^(u, 1) + C 2 {u,v)a* n (1, v)\ 

0<u,v<l 



: 0(n- 1/4 (logn) 1/2 (loglogn)) 



(5.1) 
1 

(5.2) 



and 



sup 

0<u,v<l 



a* n {u,v) 



W*(u,v) 



(uAv) 5 (l -uAv) s (uAv) s (l -uAv) s 
Let A™*(t; A) denote the solution to 



= Op(l). 



(5.3) 



{C* (u 1 "', u*) - u a }X{u, t) du = 0. 



Then (5.1) implies that 



{AT(t; A): < t < 1} i {A%(t; A): < t < 1}. 



(5.4) 



Write 
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{-u A ^}\(u,t)du 

(n/(«+l)) 1 / ( < 1 " t > vt ) 



{C^u 1 - 1 ,«*) - u A{t) }X{u,t)du 

(n+lj-i/td-tjvt) 



A(t)i 



{l-u AW }X{ Ul t)du 

(n/(n+l)) 1 /((i-t)vt) 

=:/i(i)+J 2 (t)+/ 3 (t)- 
Because 1 > A(t) > (1 - t) V f > 1/2, (2.3) implies that Ji(t) and I 3 (t) are finite 

.(n+l)- 1 /" 1 -'^') 

sup \/n\h(t)\< sup v 7 ™ / w 1 / 2 A(u,t)du = o(l), 

o<t<i o<t<i Jo 

sup Vri\l3(t)\< sup \fn \ (l — u)X(u,t)du = o(l). 

0<t<l 0<i<l J(n/(ri+l)) 1 /((i-*)vf) 



From the condition 



sup 

o<t<i Jo 



{ u ( 1 -t)v*(i _ u (i-t)vtjj& A(M) du < , 



in (2.3) and (5.3), it follows that 

rl 



sup 

0<t<l 



(g^(u , u j - W* (u L -\ u t ))A(u, t) du 



By (1.1), we have 

< C^u 1- *,***) = u A ^-^{A{t) - tA'(t)} 
< u ^-m-^-t){ A {t)-tA'{t)} 

and 

< C 2 {u 1 - t 1 u t ) = u^-^Ait) + (1 - t)A'(t)} 
< u ( 1 - t ) vt -*{A(t) + (1 - t)A'(t)}. 

Because A(t) and -A'(t) are bounded on [0,1], from the conditions 

sup / ,1 u (i-*)v*-(i-t) u (i-t)5 2 ( 1 _ „!-*)*> A(tt,t)du < oo 

0<t<l Jo 

and 

sup / u ( - 1 - t '> vt - t u tS2 (l-u t ) S2 X{u,t)du<oo 



o<t<i Jo 
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in (2.3), it follows that 



sup 

0<t<l 



sup 

0<*<1 



<(u 1_ \ l)Ci(u 1-t , u*)A(u, t) du 
1 

W*(u 1 - t ,l)C 1 (u 1 - t ,u t )X(u,t)du 
<(1, u*)C 2 (u 1 - t , u*)A(u, i) du 
W*(l,« t )C , 3 (tt 1 ~ t ,«*)A(« I t) du 



(5.8) 



By the condition 



sup n 

0<t<l 



-l/4+5i 



(«/(>l+l)) 1/((1 " t)vt) 



'(„+l)-i/((i-t)vt) 

in (2.3), (5.2), (5.7) and (5.8), wc have 



sup 

0<t<l 



X(u, i) du — > 



VEl 2 (t)- I {W*( U 1 -* ! u*)-W*(u 1 - t ,l)C 1 (u 1 -*,u*) 
-W*(l,u t )C 2 (u 1 - t ,u t )}X(u,t) du 



(5.9) 



O p ( n~ 1 ' 4 (log n) 1/2 (log log n) X ' A sup 



(n/(n+l)) 1/((1 - t)vt) 



0<t<l J (n+1) vt) 



A(u,t)du) +Op(l) = Op(l). 



By (5.6) and (5.9), we have 



sup 

0<t<l 



v /^{C'*(u 1 - t ,u t ) -u A(t) }A(u,i)du 



{W-fa 1 , u l ) - W*(u L -% IJCiCu 1- *, u J ) - u*)C2(u , u*)}A(u, t) du 



= o P (l), 
which is equivalent to 



sup 

0<t<l 



^/r(y)_/(i)} A ( U)t ) du 



(5.10) 



W*(l,u t )C 2 (w 1 - t ,u t )}A(M,^du 



= Op(l). 
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The foregoing equation shows that as n — > oo , 



sup 

0<t<l 



{ 



( ^T(t;A) 



u AW }X{u,t) du 



o p (l), 
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(5.11) 



which implies that 
P(A%*(t;X) > 4/3 for some t e [0, 1]) 



<p 



sup 

Vo<t<i 



> inf 

0<t<l 



\ u Mt) - u 4 / 3 )X{u,t)du 



since 1/2 < A(t) < 1 for all < t < 1. Similarly, 
P(A™*(t; A) < 1/3 for some t G [0, 1]) 



<p 



sup 

Vo<t<i 



> inf 

0<t<l 



(u^ 3 -u A ^)X(u,t) du 



— > as n — > oo. 
Thus, 

P(l/3 < A™*(t; A) < 4/3 for all i e [0, 1]) 1. 
By the mean value theorem, 

" -t/W}A(u,t)du 



i a( U ,t)A( t ) + (l-a(n, t ))i»*( t ;A) (logw)A(w ^ )du 



(5.12) 



(5.13) 



x(#)-4r(t;A)) 



for some a(u,i) £ [0,1]. Because 1/2 < A(t) < 1, we have < a(u,t)A(t) + (1 - 
a(u,i))A™*(i; A) < 7/3 when < i™*(t; A) < 4/3. Thus, from (5.12), it follows that 



P| inf 

. 0<*<l 



> sup 

0<*<l 



a( U ,tM(t)+(l-a(u,t))^*(t i A)(_ logtt j A(ttjf j du 

f u 7/3 (-\ogu)X(u,t)du^ -> 1 



as n—too, which, combined with (5.11), (5.13) and (5.1), implies that 

sup \Ar(t;X)-A(t)\=o p (l). 

0<t<l 



(5.14) 



Then sup 0<4<1 \A%(t;X) - A(t)\ =o p (l) follows from (5.14) and (5.4). 
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We next prove that A™(t; A) is continuous for t G [0, 1]. For t m ,t G [0, 1] and t m — > t G 
[0, 1] as m — > oo, we have 

1/2 . /-I 

u *«(*m;A)A(u,i m )du + / ( u A n (t m ;\) _ i)X(u,t m )du 
J 1/2 

1/2 ,1 

C„(u 1 ~ tm ,u tm )\{u,t m )du+ / (C„(u 1 -* m ,u tm ) - l)X{u,t m )du. 
la Ji/2 

Note that the function 

,1/2 ,1 

C n (u 1 - t ,w t )A(w,t)du+ / (C„(u 1 - t ,w t )-l)A(u,t)dw 

Jl/2 

is continuous in t G [0, 1]; thus, we have 

lim f / 7 w 1 ™(* m;A )A(u,t m )du+ / (w i "(* m;A )-l)A(M,t m )d^ 
m ^ 00 \Jo Jl/2 / 



1/2 

1/2 r l 

C n (u 1 - t ,u t )X(u,t)du+ / (C , „(u 1 ^ t ,M t )-l)A(M,t)du. 

Jl/2 



Because 



1/2 ,1 

w Q A(u,i)dw + / (u a -l)A(u,t)du 

Jl/2 

is continuous in i G [0,1] and is monotone in a for each t G [0,1], we conclude that 
4j?(t m ; A) -> A) as to -> oo. Thus, i™(t; A) is continuous in [0, 1]. 

Note that 

sup f\u^ vt (l - u^ vt )} 52 X(u, t) d« < oo 



0<t<l JO 



for some ($2 G [0, 1/2) in (2.3) implies that J^ 2 i/ 2 X(u, t) du < oo. Thus, using 

u a(u,t)X(t)+(l-a(«,t))^ , (u;A)(_l g tl ) A ( tl>t ) 

= u A ( t )(-log U )A( u ,t) w ( 1 - a ("' t »^*^ A )- A ( t » 
<^W(-logu)A(u,t)u-( 1 -°^ i » sup os*<il i n*(*'- A )- A WI, 

A(t) > 1/2 for all tG [0,1], (5.14) and 

< ir Sl — 1 < — u~ S2 for all m G [0, 1] and any fixed < si < s 2 < 1, 

S2 
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we get that 



sup 

0<t<l 



< sup 

0<t<l 



< sup 



«(«.*)A(t)+(l-a(« > t))A:«(t;A)( logu ) A(U)t)du _ / (l og u ) A ( u , t) du 

JO 

^(-log^Xiu^t)^-*- 1 -^ 11 '^^ ^ \AT^,X)-A(s)\ _ 

(1 - o(«,t)) su Po<s<1 |ir (*; A) - A(s)\ 



o<t<i V (1 - a) supo^^! |^*(s; A) - A{s)\ + (A(t) - tf 2 )/2 

x f U -( 1 -»K*))™Po<a<il^*( S ;A)-A( s )|+(A(t)+5 2 )/2 ( _ logu)A(U;t)du 



(5.15) 



= o p (l)OJ sup / u -( 1 -°( u '*)) su Po<«<il^*(^)-^WI+(^(t)+*2)/2(_iog U )A( u ,t)du 
\0<t<l Jo 



= o p (l)O p sup / w' )2 (l-u <)2 )A(u,t)dM =o p (l). 
\o<t<i7o / 

Note that the two processes A™(t; A), B(t) are continuous for t E [0, 1]. Thus, from (5.1), 
(5.10), (5.13), (5.15) and (5.4), we conclude that y/ri{A%(t; A) - A(t)} converges to B(t) 
inC*([0,l]). □ 

Before proving Theorem 3.1, we present some lemmas. Throughout, we assume that t 
is a given point in (0, 1) and use 9o to denote Ao(t). 

Lemma 5.1. Under the conditions of Theorem 3.1, as n— > oo, we have 



Proof. Write 

Vi{u,t) = K 



A 



WO-^Caiu 1 -*, u')}A(u, i) du 



u - F 1 n [ {1 ~ t \x ]1 ) \ K ( u - F^(X j2 ) 



« F^iX^ fu- F^_ t {X 32 ) 
A 



? i/t 



(5.16) 



V a (u,t) + Vi2(u,t) 
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and 

n 

n z — ^ 



n 

i=l 



n 



/l-b„ n 
Y,{Vii{u,t)-u e }\(u,t)du 
i=l 

/l-b n " 
^K; 2 (u,t)A(u,i)du 
1=1 



n" 1 



. i=l 



n- 1 



i=\ 3=1 



xK ( ^- p n?( x *) y (%t)du (5 . 17) 



EE*( 



x ^ K( M ) ^ ^ ) ^(u,f)d« 



eew wi -r 

»=i i=i 



=:Ii + ---+7 4 - 
Furthermore, the first term I\ can be expressed as 



h = A(M) / / -E / (^l,-i(^l)<Sl,^n2,-i(^2)<S2)/l- 



h 

'" fc ^'){j[ 1 j( 1 ;t<*-™^(- + ;)' 



x M "--f"''\ f H_£^') d,,!'"-" d4" - „»° } d« (5.18) 
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F n a(X a )<——[82 



n + 1 \ n 



x h~ 2 k 



,i/(i-t) 



x k 



--^\ ds {/( 1 -»ds 1 2 /t -u»°'>du 



i-6n r 1 r 1 



lJ-i 



, ,i . 1 n . , . f 1 
(u-s 2 /i) + 



A(n,i)^C„ ——(u-8 1 h) 3 - t + 



n+l 



n + 1 ?i + l 
1 



n+l 



— (u-^/i) 1 "^ 
1 n + 1 n + l 



(u-s 2 /i) + 



n + l 



n + l 



fc(si)fc(s 2 ) ds\ ds 2 du 



1— 6n /•! rl 



A(M)<j C( -4-(u - s^) 1 "* + - s 2 h) t + 

'n + l n+ln + 1 



- C{v} t ,n t )|fc(si)/c(s 2 )dsids 2 dn 

=://i + 77 2 . 
Because sup an<u<1 _ hn (ft/u) < ft/tin ~~ > an d 

inf min{(n + l)u*, (n + l)u 1_t } > (n + l)a„ -> oo 

a n <u<l—b n 



n + l 



as n — ► oo , we have 



sup 

, <«<!-&„, -Ks<l 



log(u — sh) 



logu 



2/i/u 2/i 
< sup — - — < 



2/; 



a n <«<i-6„ -logu -a„log(a n ) -(l-6„)log(l-6„) 

1 



-+0, 



and 



sup 

,<«<1-L,-Kj<1 



sup 

,<u<l-6.,-Ks<l 



n+l 



(u-sh) l - t + 



n + l 



- 1 



n+l K n+l J 



- 1 



(5.19) 



-+0 (5.20) 



0, (5.21) 
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which, together with (1.1), imply that 

log 
+ log 
A (log 



sup 

.<u<l-6„,-Ksi,s 2 <l 



(u-sxhy S + 



77+1 



sup 

a n <.u<l — b n , — 1<si,S2<1 



n+ 1 



4 


71+1 


t _ 


1 




1 - 

n+ 1 


+ 


1 1 



3gU — 1 



71+1 



lot 



71+1 



(u-sih) 1 -^ 



71+1 



77 , .. t . 1 



+ log<^ — — -(u-aaft)* + 

1 71 + 1 71+1 



sup 

a ra <u<l — 6 TlJ — 1<si,S2<1 



sup 

a n <u<l — b n , — 1<si,52<1 



sup 

a„<u<l-6„,-l<si,s 2 <l 



sup 

a™ <u<l — 6 n ,— 1<si,s 2 <1 



sup 

a n <u<l — b n ,— 1<si,S2<1 



sup 

, <U<l-fc„,-l<Sl,S2<l 



Ci 



r* 2 



Ci 



n > Ml-t 1 



71+1 

n 

71+1 

n 
n+ 1 
n 

71+ 1 

n 

71+ 1 

n 



( u -s 1 hy- t + 

(u-s 2 hf + 



71+ 1 



77+ 1 



-A(t) 



(u - s^) 1- ' + 
(u-s 2 ft)* + 



77+ 1 



1 



71+ 1 



->o, 



(u-«ifc) 1 -* + 



77+1 



^-(u- S2 h) t + ^-]-C 2 {u 1 -\u t ) 

77+1 71 + 1 / 



■ \1-f 

u-sih) 1 f + 



c 



12 



C 2 



n 


+ 


1 




n 




u 


+ 


1 




n 




ii 


+ 


1 








n 


+ 


1 








ri 


+ 


1 




n 




77 


+ 


1 



77+1' 



U - s 2 hf + 



77 + 1 



77. -Sl/i) 1 ' + 
77- S 2 /i)* + 



77+1 



77+1 



u- s 2 hf + 



77+1' 



77+1 



C^Ti 1 -*,^) 



0. 



(5.22) 



Thus, by (3.2), (5.22), and similar arguments used in the proof of Theorem 2.1, we can 
show that 



Vn//l4 / {l¥(7t 1 -^7t t )-Fy(7i 1 -^l)Cl(7t 1 -^7t t )-VK(l,7t t )C'2(7t 1 - t ,7t t )}A(7i,^)d7i. (5.23) 
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It is straightforward to verify that 



\C 2 {u l - t ,u t )u t \=0{u A W)=0{u 1 ' 2 ), 
ICi^u 1 - 4 ,^)^ 2 - 24 !! - logu}| = 0(u A(t) ) = Oiu 1 / 2 ), 
lOafa 1- *, «*)«** logtt| = 0(w A(t) ) = 0(w 1/2 ), 
I ICiaCu 1 -*, «*)«{! - log U }| = Ofa^) - O^ 1 / 2 ) 



uniformly for u 6 [a„, 1 — b n \. By Taylor's expansion, we have 
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(5.24) 



l-i>n /"l /•! 

lJ-1 



n+1 



1 - 



+ C 2 (u 1 - t ,M*)w* 



1 



n+1 



Sl/l 



s 2 /i 



l-t 



+ -C 11 (u 1 -\u t )(l + o(l))u 2 



n+1 



1 - 



Sift. 



l-t 



1 



(n+ 



(n + l)u 



(n + 



- 1 



l-t 



- 1 



+ 2^22(" 1 - t ,^)(l+0(l)) W 2t 



n+1 



s 2 /i 



1 



+ Ci 2 (w 1 - t ,u*)(l+o(l))u 



s 2 h\ t 



n+1 



n + 1 
1 



s\h 



(n + l)u* 
1 

(n + l)it 1-: 



(5.25) 

2 

1 ' 



u ) (n + l)u* 



x fc(si)fc(s 2 )A(u, t) dsi ds 2 du. 

Consider the first term in the foregoing expression. By (3.2), (5.22), (5.24), and the 
symmetry of k(s), we have 



n + 1 



s\h 



(n + l)u x -* 



x fc(si)fc(s 2 )A(u, i) dsi ds 2 dn 



1 — 6„ 



n + 1 



n + 1 



1 - 



1 \k(si)X(u,t)dsidu 



l-b„ 



Ci(u 1 -*,w t )A(u,t)dw 
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(l + o(l)) 
1 



l-6n d 



Ci^ 1 "*, u*)?/" 1 -*— -(1 - t)(-t)sik( Sl )X(u,t) dsi du 

J-i 2(n+l) 



n + 1 



l-b„ 



C 1 (u 1 -\u t )(l-u 1 - t )\(u,t)du 



Oh 



l-b„ 



3/2 X{u,t)du) +0 1/2 \{u,t)du) =o(l/V») 



Other terms of (5.25) can be handled in the same way, resulting in 



ff 2 = o(l/Vn)+0 



i-&„ 



1 



l-b„ 



l-6„ 



c 11 (« 1 - t J « t )« 2_2t (- + 



u 2 (n + l)tt 

2 



A(u,i) 



du 



1 



u (n + l)u 



X(u,t) 



du 



C 2 2(u 1 - t ,U t )u M [- + 



1 



u (n + l)u 
1 



A(it,t) 



u (n + \u (n+l)u 



1 



l-b„ 



o(l/\AI) + /i 2 / u~ 3 / 2 A(u,i)du 



A(u,t) 



+ o(h 2 J^ " {\ogu}- 1 u-^ 2 X{u,t)du 
= o(l/V^). 

For the second term, in (5.17), by the mean value theorem, we can write 



du (5.26) 



l-b„ n n f Al/(l-t) 



EE 



gxg^ - FT- t} (X 3l ) Ju- F^\X 3l ) 



(5.27) 



X i£ 



l/t 



fe)\ 



X(u, i) du, 



where £n,i,j is between F n i(Xji) and i* , n i ) _j(.X J -i). Using the equation 
F n i t -i(Xji) — F n i(Xji) = —F n i(Xji) — —I(Xii < Xji), 
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we have 

sup ihiix^) ~ F n iM x n)\< n ~\ 

l<i,j <n 



sup li^ 1 -^!) - ^, l ^(Jfii)l < T— 

Cz.7<n ^ t 



V(i-t), 



Then, uniformly for u € [a n , 1 — 6 n ], 
-i/Ci-0 
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(5.28) 



EE 7 

i=i j=i 



< i 



n n / 1 

< EE P ( ( u - /i ) 1 ~ t - - ^ ^i(^i) ^ ( w + + 



i=i j=i 



< n x 



(n + l)(u + hy-t + (n + l)/n (n + l)(u - h) 1 - 1 - (n + l)/n - 1 



(5.29) 



:0(it _1 n/i) 



and 



E J 

j=i 



1/(1-*) 



1(1 



< l ) =0(ur 1 /i). 



(5.30) 



Because k(s) is a density function with support on [—1,1], it follows from (5.27), (5.29) 
and (5.30) that 



h = 0\ h~ L n 



+ 



( h~ 2 n 



l—b n n n 

EE 7 

i=i j=i 

l-b„ " n 

EE 7 

i=i j=i 



u-F. 



V(i-t) 



POO 



< 1 )A(u,*)du 



< 1 )X(u,t)du 



(5.31) 



0( n" 1 



i-fe„ 



1 A(u,t)du I = o(l/Vn). 



Similarly, we can show that 

7 3 =o(l/\/n) and I4 = o(l/ y/n) . 
Thus, the lemma follows from (5.23), (5.26), (5.31) and (5.32). 

Lemma 5.2. Under conditions of Theorem 3.1, we have 

1 " / r 1 

-E^°) 4i n / Wtt 1 -*,"*)-^ 1 "*.!)^^ 1 -*,^) 



(5.32) 

□ 
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_ 2 

- W(l, u^C^m 1 "*, u*)}A(u, t) du ' 

as n — > oo . 

Proof. By (5.16), we can write 

QKo) 

{Vn (ux, t)V a («a , t) + Vji («i, t) V i2 («a , «) - $i («i , «)«a 

+ V a («i , t) Vii (« a , t) + Vi 2 («i, t) V l2 (u 2 , t) - V a (ui , t)u 9 2 
- UiVa(u2,t) - ulV i2 (u2,t) + u e 1 U2}X(ui,t)X(u2,t)duidu2- 

Using arguments similar to those in (5.27), we have 

n /-l-bn rl—h 



i=l 



/ / V r i2(wi,0^2(w2,i)A(ui,t)A(u 2 ,i)dMidu 2 

1 V^V^ f Fnl,-i{Xj\) -F nl (Xji) 1 £,t/(l-t) fv 



a n J a 



n 

K i=i j=i i=i 



h , u l -F^ 1 - t \x ]1 ) \ K fu l -FlJ 2 t {X ] 2) 



h J \ h 

Fn2,-i(Xj2) - F n2 (X j2 ) 1 A(l-t)/t/ y x 
ft ~t F ^ {X]2) 

^-F 1 n l2 t {X ] 2) \ K ( u 1 -F 1 n [ {1 - t) ^ 



h \ h 



F„i -i(Xn) - F nl (Xn) 1 #,i/(i-t), v \ 

h Y^~t nl {Xll) 

• U2-F 1 n ^- t \x ll ) \ K (u2-F 1 n i\X l2 ) 



(5.33) 



h 

x A(ui, t)X(u2,t) du\ du2 + o p (l). 



h ) \ h 

F n 2,-i(Xl2) -F n 2{Xl2) 1 »(i_ t )/ t 
+ £ ^„2 
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It is straightforward to check that 
1 " 

-J2{Fm(x) - F nl! _i(x)}{F nl (y) - F„i,_,(i/)} 



n 
»=i 



n + 1 » , . n + 2- - 
— =-F„i(a;Ay) ^-F nl (x)F nl (y), 

1 " 

- J^{F n2 (x) - F nai _ < (a:)}{F„2(i/) - F n2 ^ t (y)} 



n . 
?— l 



and 



n + 1 - , . n + 2 - , . - . . 
— Ay o-F n2 (x)F n2 (y) 



1 n 

- ^{F„i(x) - F„i,-,(a;)}{F„ 2 (y) - F„ 2 ,_;(y)} 



n 

i=l 



n 

Then (5.33) can be written as 

n r l-b n rl—b. 



\c n {F nl {x),F n2 (y)) - ^^F nl {x)F n2 { y ). 



^ ,L pi — On r l — n 

/ / Vi 2 (wi,t)Vi2(u2,*)A(ui,t)A(u2,t)duidu 2 

n i=1 Ja„ Ja„ 



J a j a [^EE{( # «i(^ AX 'i)--p»i(^i)^i(^x))- 



/ £nl/(l-t)/l 



x A' 



(42(^-2 AI, 2 )- F n2 (X, 2 )F n2 (X l2 ))^ 



x X 
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+ {C n (F n i(Xji),F n 2(Xi2)) — F nl (Xji)F n 2(Xi2)) 

b n\ \ X 3l) t n2 \ A W 



t(l-t) 



U2-F^ t (X l2 )\ T Ju2-F^ 1 - t \x ll ] 



x jfc — — r v ' \k 



+ (C n (F n i(Xii), F n2 {Xj2)) — F n i(Xii)F n 2(Xj2)) 

pl/(l-t)/ y \pM/Vy \ 



*(!-*) 
x k 



u 2 F^ 1 ~ t) (X ll )\ r Ju 2 - F^(X l2 ) 



x fc 



ft 

%^) V( ^-& /(1 -°( 

ft J \ h 

x A(iti,i)A(u 2 ,i) dui du 2 + o p (l). 

Based on the foregoing decomposition, we can show that 

fl-6„ /■!-&„ 



- ^ J / / V l 2(u 1 ,t)V l 2(u2,t)X(u 1 ,t)X(u 2 ,t)du 1 du 2 

j=l J^n Ja n 

L 

({«i _t a 4"* - ul~*u5~*}Ci i4)ai(t4 -t ,t4) 
+{c(«i- i ,4)- u ;-^}c 1 («|- t , !1 i)c 2 K- ! 1 ^) 

+ {C(t4-*, ui) - t4-*ui}Ci (u 2 - 4 , w 2 )C 2 , ui)) 
x A(w,i, i)A(u 2 ,i) dui du 2 + o p (l). 



(5.34) 



Similarly, we have 

rl-b„ r l-b n 



- ^ / / V r a(wi,0^2(u 2 ,i)duidit 2 



{c(4~*, 4) - c(«i _ * a u\-\ u[)d (t4~*, 4) 

+c(4-^u t 1 ) U2 c' 2 ( U2 -^ U2 )-c'(4-^u^A U2 )c 2 ( U2 -^ U2 )} 

x A(ui,t)A(u2,t)duidu2 +o p (l), 



(5.35) 
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n rl-b n t-l-br, 



/ / V il (u l7 t)V ll (u2 1 t)\(u 1 ,t)\(u2,t)du 1 du 2 



1 r l 

C(wJ _t Au^.iii Aul)X(ui,t)X(u2,t)duidu2 + o p (l), 
- ^ / / V r a('"i,<) u 2 ^( M i^)^('"2,i)duidu2 

C(wi"*,M' 1 )C(u2" t ,i4)A(iti,t)A('U2,i)dwidM2 +o p (l), 



(5.36) 



(5.37) 



and 



~Y. / V r 4 2(ui,0«2 O A(wi,i)A(w 2 ,t)dw 1 dw 2 -:o p (l). (5.38) 

Thus, the lemma follows from (5.34)-(5.38) and the fact that 

e( f {Wiu 1 ' 1 ^ 1 ) - Wiu 1 -*, l)Ci(« 1_ *,u*) - ^(l,w t )C* 2 (u 1 ~ t ,M t )}A(u,i)dw 



1 r l 

{C(wJ-* A t4~*i«i A u\) - C^'VlM^*^) 

- (C(«|- f A ^'X) - C^-Vi)^"')^^-^) 
-(Ciul-^ulAuD-Ciul-^uDul^ul-^ui) 

-{Ciu^ Aul-\u\- t )-u\- t C{u\-\u2))C 1 {u\- t A) 

+ (c( ll ;- t I 4)-^- t 4)c 1 (ui- < ,« ( 1 )c 2 ( l 4- t ,4) 

- (C(*4~*,Ui Aw 2 )-w t 1 C(^-*,w*))C*2(ul" t ,u' 1 ) 

+ (c(^- t , u t 1 )- W2 -V 1 )C2K- t , w t 1 )c 1 (^- t , M2 ) 

+ (u'j Au 2 - u\ul)C2(ul~\u\)C2(ul~\u 2 )}\(ui,t)\(u2,t) duidu 2 . □ 



o Jo 



Proof of Theorem 3.1. Using similar expansions as in the proof of Lemma 5.1, we can 
show that maxi<i<n |Qj(^o)| = (^(n 1 / 2 ). Thus, using Lemmas 5.1 and 5.2 and standard 
arguments in expanding the empirical likelihood ratio (see, e.g., Owen [13]), we obtain 
that as n — > oo , 



28 



L. Peng, L. Qian and J. Yang 



Acknowledgements 

We thank the editor, an associate editor and two reviewers for their helpful com- 
ments. Peng's research was supported by NSA Grant H98230-10-1-0170 and NSF 
Grant DMS-10-05336. Qian's research was partly supported by the National Natu- 
ral Science Foundation of China (Grant 10971068). Yang's research was partly sup- 
ported by the National Basic Research Program (973 Program) of China (Grant 
2007CB814905) and the Key Program of National Natural Science Foundation of China 
(Grant 11131002). 

References 

[1] Bucher, A. and Dette, H. (2010). A note on bootstrap approximations for the empirical 
copula process. Statist. Probab. Lett. 80 1925-1932. MR2734261 

[2] Bucher, A., Dette, H. and Volgushev, S. (2011). New estimators of the Pickands 
dependence function and a test for extreme-value dependence. Ann. Statist. 39 1963- 
2006. 

[3] Caperaa, P., Fougeres, A.L. and Genest, C. (1997). A nonparametric estimation pro- 
cedure for bivariate extreme value copulas. Biometrika 84 567-577. MR1603985 
[4] Chen, ,]., Peng, L. and Zhao, Y. (2009). Empirical likelihood based confidence intervals 

for copulas. J. Multivariate Anal. 100 137-151. MR2460483 
[5] Deheuvels, P. (1991). On the limiting behavior of the Pickands estimator for bivariate 

extreme-value distributions. Statist. Probab. Lett. 12 429-439. MR1142097 
[6] Falk, M. and Reiss, R.D. (2005). On Pickands coordinates in arbitrary dimensions. 

J. Multivariate Anal. 92 426-453. MR2107885 
[7] Fermanian, J.D., Radulovic, D. and Wegkamp, M. (2004). Weak convergence of em- 
pirical copula processes. Bernoulli 10 847-860. MR2093613 
[8] Genest, C. and Segers, J. (2009). Rank-based inference for bivariate extreme-value cop- 
ulas. Ann. Statist. 37 2990-3022. MR2541453 
[9] Gong, Y., Peng, L. and Ql, Y. (2010). Smoothed jackknife empirical likelihood method 
for ROC curve. J. Multivariate Anal. 101 1520-1531. MR2609511 
[10] Hall, P. and Tajvidi, N. (2000). Distribution and dependence-function estimation for 

bivariate extreme-value distributions. Bernoulli 6 835-844. MR1791904 
[11] Jing, BY., Yuan, J. and Zhou, W. (2009). Jackknife empirical likelihood. J. Amer. 

Statist. Assoc. 104 1224-1232. MR2562010 
[12] Kojadinovic, I. and Yan, J. (2010). Nonparametric rank-based tests of bivariate extreme- 
value dependence. J. Multivariate Anal. 101 2234-2249. MR2671214 
[13] Owen, A.B. (1988). Empirical likelihood ratio confidence intervals for a single functional. 

Biometrika 75 237-249. MR0946049 
[14] Owen, A. (1990). Empirical likelihood ratio confidence regions. Ann. Statist. 18 90-120. 
MR1041387 

[15] Owen, A. (2001). Empirical Likelihood. New York: Chapman & Hall/CRC. 

[16] Pickands, J. Ill (1981). Multivariate extreme value distributions. Bull. Inst. Internat. 

Statist. 49 859-878. MR0820979 
[17] Qin, J. and Lawless, J. (1994). Empirical likelihood and general estimating equations. 

Ann. Statist. 22 300-325. MR1272085 



Dependence function for an extreme-value distribution 



29 



[18] Segers, J. (2012). Asymptotics of empirical copula processes under nonrestrictive smooth- 
ness assumptions. Bernoulli 18 764-782. 

[19] Zhou, M. emplik: Empirical likelihood ratio for censored/truncated data. R package version 
0.9-3-1. Available at http://www.ms.uky.edu/~mai/splus/library/emplik/. 

Received October 2010 and revised September 2011 



